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The common variants in lysyl oxidase-like 1 gene (LOXLl) axe associated with exfoliation glaucoma (XFG) 
patients developed through exfoliation syndrome (XFS). However, the risk allele of a variant in LOXLl has 
been found to be inverted between Asian and Caucasian populations. Therefore, we newly performed a 
genome-wide association study using 201 XFS/XFG and 697 controls in Japanese, and identified 34 
genome-wide significant single-nudeotide polymorphisms (SNPs) distributing in not only LOXLl but also 
TBC1D21 and PML at the 15q24.1 locus. These SNPs were confirmed by an independent population 
consisted of 121 XFS/XFG and 263 controls in Japanese. Moreover, further analyses revealed a unique 
haplotype structure only from the combination of TBC1D21 and LOXLl variants showing a high XFS/XFG 
susceptibility specific for the Asian population. Although there still should be other gene(s) in the other 
region(s) contributing to the disease process, these results suggested that the combination of newly 
discovered variants in these genes might be useful for precise XFG risk assessment, as well as for elucidating 
the molecular mechanism of XFG pathogenesis through XFS. 

Glaucoma is one of the leading causes of irreversible blindness worldv\dde'. It is a complex age-related 
disorder consisting of heterogeneous subtypes that share common features for the development of the 
disease, i.e., progressive loss of retinal ganglion cells and optic nerve axons that result in visual field 
defects^. Therefore, it is important to investigate each distinct type of glaucoma in order to elucidate the different 
molecular pathways that lead to the common optic-nerve degeneration. 

Exfoliation glaucoma (XFG) is categorized as secondary open-angle glaucoma and reportedly develops from 
exfoliation syndrome (XFS), as manifested by an observation of abnormal fibrillar deposits (known as pseudoex- 
foliation (PEX)) on the lens and iris epithelium^. In general, XFG has a worse prognosis than primary open-angle 
glaucoma (POAG), the major type of glaucoma, as it is often observed in patients as elevated intraocular pressure 
(lOP) combined with severe damage of visual function at the time of initial presentation due to the asymptomatic 
nature of glaucoma. Moreover, according to the findings reported in the Tajima Study'', a robust epidemiology 
study conducted in Japan, XFG accounts for 0.8% of the 5.0% of all types of glaucoma observed in Japanese people 
over 40 years of age, and the prevalence of XFG increases with advancing age. Consequently, it would be of 
enormous benefit if the predisposing genetic factors of XFG could be used to predict the risk of developing the 
disease prior to onset in order to stop, or at least hinder, the irreversible progression toward blindness by initiating 
early medical intervention. 

In 2007, Thorleifsson et al. reported a genome-wide association study (GWAS) using the populations from 
Iceland and Sweden, and described that one intronic variant and two nonsynonymous coding variants of lysyl 
oxidase-like 1 gene {LOXLl) confer the risk of developing XFG, possibly through XFS^. From the functional point 
of view, the involvement of LOXLl in XFG pathogenesis is quite rational, because LOXLl plays an indispensable 
role in elastic fiber homeostasis and maintaining the structures of the trabecular meshwork, the tissue that 
determines the lOP level in the eye''. However, LOXLl, in and of itself, fails to sufficiently explain the complex 
mechanism of the disease progression of XFG through XFS. This observation is also supported from the point of 
view of genetics, as many replication studies using populations derived from different genetic backgrounds 
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(Caucasians from the US' ' and Europe'"'", Asians from Japan'^""" 
and China" '", and Africans from South Africa") have succeeded in 
replicating the association of the LOXLl variants, although the allele 
frequency of the variants varied among the ethnicities. In fact, the 
risk allele of one of the exonic variants (rs 1048661) has been found to 
be inverted between Japanese'^ and Nordic populations^. Conse- 
quently, it is wise to investigate other modifying or causative variants 
in genes other than LOXLl by performing a GWAS using a non- 
Caucasian population in order to identify variants that are determin- 
ing the different etiology of XFS onset and its progression to XFG 
among the ethnicities. 

In this present study, we conducted a GWAS and replication study 
using two independent populations of Japanese case-control subjects 
in an attempt to identify new variants associated with XFS/XFG. As a 
result, we identified a cluster of genome-wide significant SNPs only 
from the chromosome 15q24.1, which is well known as a locus 
including LOXLP. However, the genome-wide significant SNPs were 
distributed in not only LOXLl, but also in the two adjacent genes, 
TBC1D21 and PML, without being in significant linkage disequilib- 
rium (LD). We also identified a suggestive association only between 
TBC1D21 and LOXLl variants by a conditional analysis, and 
revealed a unique haplotype structure of the variants derived from 
these two genes in our study population. Together with the accu- 
mulating knowledge of LOXLi -XFS/XFG relationships, the variants 
and genes newly discovered in this study would help to provide a 
better prediction of XFG risk assessment, as well as to elucidate the 
molecular mechanism of XFG pathogenesis through XFS. 

Results 

GWAS for XFS/XFG. The samples used in this study are sum- 
marized in Table 1. In order to perform GWAS, we first selected 
the XFS/XFG subjects based on our strict diagnosis criteria, and 
then obtained the genotype data. We performed a GWAS using a 
population consisting of 201 XFS/XFG patients and 697 healthy 
controls (Table 1). In total, 652,792 SNPs that passed the quality 
controls (QC) were used for the association study (Figure lA). 
According to the quantile-quantUe (Q-Q) plot (Supplementary 
Material, Figure SI), the genomic inflation factor (/t) showed 1.020, 
suggesting that the population substructure should not have any 
substantial effects on the association analysis. Under those 
conditions, we observed a strong association signal from a locus on 
chromosome 15 (Figure 2A). In fact, we obtained 34 genome- wide 
significant SNPs (Supplementary Material, Table SI) that passed the 
Bonferroni correction threshold (0.05/652,792 = 7.66 X 10"") on 
chromosome 15q24.1 (Figure 2A). Surprisingly, these significant 
SNPs were found to be broadly distributed in not only LOXLl, but 
also the adjacent two genes, TBC1D21 and PML. The significance for 
this region was also confirmed by an imputation analysis (Figure 2B). 
According to the LD plots (Supplementary Material, Figure S2A), the 



Table 1 | Sample Information 



A. GWAS 



Study 


Group 


Subgroup 


n 


Age 


F/M 


GWAS 


XFS/XFG- 


Total 


201 


74.3 ± 8.9 


0.9 






XFG 


162 


75.3 ± 8.5 


0.8 






XFG with other 


21 


72.3 ± 10.5 


0.9 






glaucoma 












XFS 


18 


66.8 ± 6.6 


3.5 




Control 


Total 


697 


58.5 ± 13.4 


1.9 


Replication 


XFS/XFG'' 


Total 


121 


74.7 ± 8.6 


1.0 






XFG 


103 


75.8 ± 7.9 


0.8 






XFG with other 


18 


68.4 ± 9.8 


2.6 






glaucoma 










Control 


Total 


263 


48.6 ± 14.8 


1.5 



Affy SNP 6.0 


: 906,600 SNPs 


XFS/XFG 201 


VS 


Control 697 


Filtered : 652,792 SNPs 



■ Passed ttie Bonferroni's 
I correction threshold 
r (P< 7.66 X 10-8 ) 



34 SNPs in 15q24.1 locus 


[TBC1D21 


L0XL1 


PML ] 



B. Replication 



C. Meta-Analysis 



Sequenom MassARRAY : 43 SNPs 



XFS/XFG 121 VS Control 263 



Combined by Mantel-Haenszel test 



jConditional and haplotype analyses 



Susceptible haplotype between 
and 



TBC1D21 



L0XL1 



"Ninety XFG and 2 XFG with other glaucoma samples are overlapped with our previous study'' 
''One XFG sample is overlapped with our previous study'^. 



Figure 1 | Study design. (A) First, a GWAS was performed and 652,792 
SNPs were analyzed for 201 XFS/XFG patients and 697 controls. As a result, 
34 genome-wide significant SNPs were obtained from the genes of 
TBC1D21, LOXLl, and PML atl5q24.1 locus. (B) Next, the GWAS results 
were replicated by analyzing 43 SNPs from the 15q24.1 locus by using an 
independent population of 121 XFS/XFG patients and 263 controls. (C) 
Finally, a meta-analysis was performed by combining the two data sets by 
use of the Mantel-Haenszel test. In addition, the pairs of significant SNPs 
were assessed by conditional and haplotype analyses. The result suggested 
the susceptible haplotype might exist for XFS/XFG from TBC1D21 and 
LOXLl variants. 

SNPs in each gene seemed to be derived from 3 distinct LD blocks, 
which was also confirmed by the LD values using 34 genome-wide 
significant SNPs (Supplementary Material, Table S2). 

Replication study. Next, we attempted to replicate the GWAS results 
by analyzing an independent Japanese population consisting of 121 
XFS/XFG patients and 263 controls (Table 1). We assessed a total of 
43 SNPs at 15q24.1 locus, including 34 genome-wide significant 
SNPs from our GWAS results and 2 nonsynonymous coding SNPs 
(cSNPs) (rsl048661 and rs3825942) in exon 1 of LOXLl that 
reportedly are strongly associated with XFS/XFG^. As a result, 
almost every genome-wide significant SNP, except 1 SNP 
(rs723434), passed the Bonferroni correction threshold in the 
replication study (0.05/43 = 1.16 X 10"^) (Supplementary 
Material, Figure S2B and Table S3), confirming the association 
results obtained by GWAS. We also confirmed the association 
result of the reported 2 cSNPs in LOXLl using our Japanese 
population, of which the risk allele of rs 1048661 was indeed 
inverted from the reported result with Caucasian. When the 
combined P values of GWAS and the replication study were 
calculated by use of the Mantel-Haenszel test^", the level of 
significance of almost every SNP increased. In fact, rs893818, 
which was the most significant SNP in intron 1 of LOXLl, reached 
to P = 8.21 X 10"** (Supplementary Material, Figures S2C and S2D, 
Table S4). The rs893818 SNP was clearly independent from the SNPs 
in TBC1D21 and PML as in reference to the LD values calculated by 
applying the whole samples used in this study (Supplementary 
Material, Table S2). Moreover, the LD blocks based on the whole 
samples demonstrated that LOXLl, TBC1D21, and PML were in 3 
distinct LD blocks, suggesting that the variants from different genes 
were independently contributing to the disease (Supplementary 
Material, Figures S2E and S2F). 
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Figure 2 | Association results. (A) The SNPs with a strong association signal from the GWAS result appeared to exist as a cluster on chromosome 1 5 (red 
arrow). Horizontal dotted line represents Bonferroni correction threshold (P < 7.66 X 10"*). (B) Because the genome-wide significant SNPs were 
intensively identified on the 1 5q24. 1 locus, we also imputed the SNPs across the locus and plotted against the positions based on the NCBI Build 36 
coordinates. The LD value of is referred to the most significant SNP (rs893818) in the GWAS. The cSNP (rsl6958445) and two genes {TBC1D21 and 
LOXLl) used in the conditional and haplotype association analyses are highlighted by red. The genetic recombination rates (cM/Mb) estimated by 
HapMap Project" in Release 22 are indicated by the thin blue line. 



Meta-analysis of two studies. We then performed a conditional 
analysis for 34 genome-wide significant SNPs on the basis of 
rs893818 according to our GWAS result. Overall, it turned out that 
the effect of rs893818 was very strong, while the variants in TBC1D21 
and PML were generally weak, although we found a suggestive signal 
(P = 0.027) from a nonsynonymous variant in exon 4 of TBC1D21 
(rsl6958445) (Supplementary Material, Table S5). These two SNPs 
were considered to be truly associated with XFS/XFG because there 
were no confounding factors such as age and gender (Supplementary 
Material, Table S6). Consequently, we obtained 3 susceptible cSNPs, 
including 2 reported cSNPs in LOXLl, showing a weak correlation 
(r^ < 0.2) to each other (Supplementary Material, Figures S3). 



Therefore, we then evaluated the combinational effect of 3 cSNPs 
(rsl6958445 from TBC1D21, and rsl048661 and rs3825942 from 
LOXLl) in order to analyze the susceptibility to XFS/XFG by 
performing a haplotypic association analysis with the samples used 
in the replication study (Table 2). As a result, haplotype H6 
(rsl6958445: G, rsl048661: T, and rs3825942: G) was most 
significant for the association with XFS/XFG. Interestingly, 
haplotype H6 was more significant than haplotypes H3 and H4, 
including the significant haplotype (rsl048661: G and rs3825942: 
G) in the LOXLl for the Caucasian population^. In addition, 
according to the dbSNP database, rsl6958445 in TBC1D21 was 
found to be monomorphic for the G allele in Caucasian and 
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Table 2 | Results of Haplotype Association Analysis 
Haplotypes SNP alleles" Freq. in XFS/XFG 


Freq. in controls 


Chi-square 


P-values'' 




OMNIBUS 


NA 


kl A 
INA 


1 00 A 


O R~7 \/ 1 n-25 
V.J/ X 1 u 


HI 


GGA 


0.0039 


0.1294 


31.9 


1.62 X 10-^ 


H2 


OTA 


0.0043 


0.0407 


7.7 


0.0055 


H3 


AGO 


0.0176 


0.0835 


12.2 


0.0005 


H4 


GGG 


0.0125 


0.1517 


33.4 


7.55 X 10-' 


H5 


ATG 


0.0319 


0.0675 


3.9 


0.0471 


H6 


GTG 


0.9297 


0.5272 


1 18.5 


1.33 X 10-2^ 


"Haplotypes were constructed with the SNPs in the 
''Pvalues were calculated by use of the 'hap-assoc 


following order: rsl6958445, rs1048661, 
' option of PLINK. 


and rs3825942. 







African populations. Consequently, we obtained a unique haplotype 
structure in this present study showing a high XFS/XFG suscepti- 
bility specific for the Asian population. 

Taken together, the newly identified genes in 15q24.1 locus, 
together with the unique haplotype of TBC1D21 and LOXLl var- 
iants, could add some information to explain the molecular mech- 
anism of XFS/XFG pathogenesis, which has been considered difficult 
to explain by LOXLl alone. 

Discussion 

In this study, we successfully obtained the genetic markers strongly 
associated with XFS/XFG patients by analyzing two independent 
Japanese populations totaling 1,282 subjects from GWAS and its 
replication study. As a result, we discovered the genome-wide sig- 
nificant variants in not only LOXLl, but also 2 other genes, TBC1D21 
and PML, both of which are located on chromosome 15q24.1, which 
is well known as the "LOXLl locus"^. We also found a suggestive 
combinational effect of the variants related to the XFS/XFG suscept- 
ibility only between LOXLl and TBC1D21, which seemed to be spe- 
cific to Asian population. 

In the first reported GWAS for XFS/XFG, Thorleifsson et al. 
detected a strong association signal only in LOXLl at the 15q24.1 
locus from a Nordic population^. Their GWAS was carried out based 
on the genotype data generated by using the Infinium 
HumanHapSOO BeadChip (lUumina, Inc., San Diego, California, 
USA) array, whereas we used the Genome- Wide Human SNP 
Array 6.0 (SNP 6.0; Affymetrix, Inc., Santa Clara, California, USA), 
in which the SNP coverage is more than three-fold. We precisely 
compared the probe disposition derived from each microarray as 
well as the association results from the two studies across the locus 
(Supplementary Material, Figure S4A and B). Although the number 
of QC-passed SNPs in the locus was different in the two studies 
depending on which microarrays and sample ethnicities were used, 
we were able to directly compare the results for some SNPs 
(Supplementary Material, Figure S4A; connected with green lines). 
As a result, a genome-wide significant SNP in LOXLl was repro- 
ducible in both studies. However, the genome-wide significant SNPs 
in other genes at the 15q24.1 locus, especially in TBC1D21 
(Supplementary Figure S4B), obtained in our study using Japanese 
samples were not significant in the first reported GWAS using the 
Caucasian samples. Moreover, we found that 1 of the significant 
SNPs in TBC1D21, rsl6958445, which together with 2 cSNPs in 
LOXLl seemed to be contributing to the disease according to the 
results of the conditional analysis, was a monomorphic SNP in all of 
the HapMap populations other than Asian (see dbSNP in NCBI). 

With regard to rsl6958445, the non-risk allele "A" was a minor 
allele, whose frequency was about 0.144-0.157 or 0.166 in Japanese, 
according to the dbSNP or our control samples, respectively 
(Supplementary Material, Table S4). This meant that the contri- 
bution of rsl6958445 to the XFS/XFG susceptibility might be limited 
for explaining the different susceptibility among the ethnicities. 
However, the results of the haplotype association analysis by 3 



cSNPs (H6: GTG in Table 2) demonstrated that the significance 
was higher than that by use of only 2 cSNPs in LOXLl 
(Supplementary Material, H4: TG in Table S7). In addition, the 
non-risk allele of rsl6958445 decreased the significance of the major 
susceptible haplotype in Caucasians (Supplementary Material, H3: 
GG in Table S7). Taken all together, we concluded that the discrep- 
ancy of the results from the two studies was due to the different 
genetic backgrounds of the samples, and not due to practical reasons 
such as the microarray density. 

The association results of LOXLl variants found in the Nordic 
population has been extensively characterized in many other popula- 
tions with different ethnic backgrounds. In particular, out of the three 
variants reported, two coding variants of rsl048661 (Argl41Leu) and 
rs3825942 (Glyl53Asp) in exon 1 were mostly analyzed by the rep- 
lication analysis'" (Supplementary Material, Figure S5). Although 
both of the variants were found to be associated with the disease 
broadly beyond the ethnicities, the risk alleles of rs 1048661 and 
rs3825942 were, intriguingly, inverted in Japanese's/Chinese'" and 
South African''-" populations, respectively. The inverted genotypes 
of the risk alleles among different populations may be reflecting the 
different patterns of historical recombination between these SNPs 
and the unidentified causative variants. However, the one thing that 
is apparent is that the particular amino-acid changes in LOXLl pro- 
tein were not functionally sufficient enough to be the causative factor 
for the disease pathogenesis. Consequently, it would have been wise 
to investigate other modifying or causative variants in genes other 
than LOXLl. Thus, we decided to perform a GWAS using our 
Japanese population in order to identify unique variants contributing 
to the onset of XFS and its progression to XFG. 

To date, only a limited number of disease-associated variants 
identified by GWAS were applicable to the actual diagnosis in the 
clinic setting. Thus far, a case in point is the diagnosis that has been 
applied to predict the risk of age-related macular degeneration 
(AMD)SS'S^. In order to predict the prevalence and incidence of 
AMD, multiple risk models have been developed by incorporating 
both the genetic and environmental factors, resulting in a marked 
increase of the accuracy to discriminate risk and non-risk indivi- 
duals^'*. However, in the case of XFS/XFG, the situation is completely 
different from that of AMD, due to the limited number of identified 
genetic factors as well as their characteristics. Although the preval- 
ence oi LOXLl variants in affected patients is strikingly high (>80%), 
which means that the sensitivity is high, the prevalence of unaffected 
controls is also high enough (up to 88%) to render the specificity 
being low, indicating that the simple allelic testing of LOXLl variants 
themselves cannot exclude the individuals who would not develop 
XFS/XFG^'. Although the variants identified in this study seemed to 
possess a weaker effect on the disease than the LOXLl variants have, 
we also discovered a unique haplotype specific for the Japanese 
population (Table 2). Combining the information of each allele 
and the haplotypes would help in developing a better diagnosis to 
predict XFS/XFG risk by the genetic testing of these variants. 
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Since XFS is characterized by clinically visible abnormal fibrillar 
deposits (i.e., PEX), progressive accumulation of those deposits have 
been considered to be a trigger of developing glaucoma by increasing 
the outflow resistance of the trabecular meshwork, resulting in the 
continuous elevation of lOP^. The LOXLl protein, which catalyzes 
the formation of elastin fibers and contributes to the trabecular 
meshwork structures, is thus probably involved in the XFS/XFG 
pathogenesis'". However, the exonic variants in LOXLl fail to suffi- 
ciently explain the complex mechanism of the disease (e.g., as if they 
were affecting the quality (amino-acid changes) and/or quantity 
(expression level) of LOXLl protein) due to the differences in allelic 
distribution among different genetic backgrounds and the high pre- 
valence of the variants in the healthy population^'. 

According to our GWAS result, the significant variants for XFS/ 
XFG other than the LOXLl variants were newly discovered from 
TBC1D21 and PML genes (Figure 2B, Supplementary Material, 
Figure S2). In fact, the TBC1D21 gene, which belongs to the TBCl 
domain family, was suggested to be related to LOXLl by haplotype 
association analysis (Table 2). The TBC (Tre-2/Bub2/Cdcl6) 
domain consists of approximately 200 amino acids and is considered 
to function as a GTPase activator of Rab proteins. However, the 
functional details of more than 40 members of the protein families 
in humans and mice, including TBC1D21, have to date not been well 
investigated^^. Although there are few functional analyses of the gene, 
our RT-PCR analysis revealed that TBC1D21 is expressed in the 
human retina (Supplementary Material, Figure S6). When we com- 
pared the sequence conservation of the variant (rsl6958445) among 
several species, we found that the variant was conserved within pri- 
mates (Supplementary Material, Figure S7). We also performed in 
silico prediction in order to evaluate the influence of amino acid 
change in rsl6958445 by means of SIFT^"" and PolyPhen-2''^, and 
programs predicted as "TOLERATED" and "benign", respectively. 
The results suggested that the amino acid change in rsl6958445 does 
not affect the function of TBC1D21 protein. Consequently, further 
functional investigation is necessary to elucidate whether or not 
TBC1D21 involves in the disease etiology of XFS/XFG. 

On the other hand, PML is functionally well characterized and 
would be a good target for analyzing the relationship to XFS/XFG 
etiology, although the contribution to XFS/XFG susceptibility was 
considered to be marginal than TBC1D21 and LOXLl based on our 
analysis. PML is a member of the tripartite motif (TRIM) family and 
famous for being involved in generating fusion protein as a result of 
chromosomal translocation with the retinoic acid receptor alpha 
gene, which is often associated with acute promyelocytic leukemia^**. 
The original function of PML protein is to form a nuclear body, 
which carries several components, such as pRB, p53, CBP, elf-4, 
Daxx, SUMO-1, etc^". Therefore, if the variants in PML cause inap- 
propriate formation of the nuclear body, broad biological functions 
would be affected, including transcriptional regulation, initiation of 
translation, and sumoylation, which could associate with the elastic 
microfibrUlopathy. In fact, our RT-PCR analysis shows that PML is 
expressed in the human retina (Supplementary Material, Figure S6). 
Consequently, we consider that it would be interesting to further 
investigate the involvement of PML in the XFS/XFG pathogenesis. 

Although we succeeded in obtaining some new variants and genes 
associated with XFS/XFG, it should be noted that there must be other 
gene(s) in other region(s) contributing to the disease in order to 
explain the complex mechanism of XFS/XFG pathogenesis. 
Moreover, under the current GWAS condition, we are unable to 
verify whether or not the variants were specifically related to the 
mechanism for the progression to glaucoma from XFS. Since there 
are patients who possess PEX, but never develop glaucoma, other 
factors must be involved in determining whether or not XFS patients 
develop glaucoma. Together with the variants/genes identified in this 
study and the accumulated knowledge of LOXLl, the identification 
of the remaining variants as well as the variants specific for XFS 



patients without glaucoma could provide a complete set of genetic 
factors associated with XFS and XFG, which should be useful for their 
diagnostic tools as well as for revealing the molecular mechanism of 
their etiology. 

Methods 

Subjects. All of the procedures in this study were performed in accordance with the 
tenets set forth in the Declaration of Helsinki and were approved by the Institutional 
Review Board of Kyoto Prefectural University of Medicine, Kyoto, Japan. Written 
informed consent was obtained from all XFS/XFG participants after receiving a 
detailed and thorough explanation of the study. Those participants were interviewed 
to determine their famUial history of glaucoma and the medical histories of other 
ocular or general diseases. The Japanese volunteers were recruited between March 
2005 and August 2012 at the University Hospital of Kyoto Prefectural University of 
Medicine in order to provide peripheral blood samples. The genotyping data of this 
current study for XFS/XFG patients and healthy volunteers was simultaneously 
obtained with the data of our previous study for POAG"*". From the complete data set, 
XFS/XFG cases and controls suitable for this study were selectively chosen from the 
complete data set based on the following strict diagnosis: 1) all XFG patients were 
diagnosed by slit-lamp examination for the existence of exfoliation materials on the 
pupil and surface anterior lens capsule with open angle, 2) lOP higher than 
21 mmHg, and 3) visual field defect according to the optic nerve appearance. The 
category "XFG with other glaucoma" group includes normal tension glaucoma or 
primary angle closure glaucoma with XFS. All of the patients and healthy volunteers 
were diagnosed by 3 ophthalmologists (Y.I., M.U., and K.M.) from a single institution. 
The demographics (i.e., subject age and male-to-female ratio) of all of the subjects 
used in this study are shown in Table 1. To examine the possible confounding effects 
of the age and gender of the subjects, correlations between the case and control 
samples were assessed by use of the Student's t test or chi-square test (Table 1). 

SNP genotyping and quality control for the GWAS. As described previously^", 
906,600 SNPs were genotyped for 2,126 Japanese subjects including POAG patients, 
XFS/XFG patients, and healthy volunteers by Genome- Wide Human SNP Array 6.0 
(Affymetrk), and strict QC was applied. For the GWAS, 203 XFS/XFG patients and 
718 controls were initially selected from the total number of subjects based on the 
latest diagnosis. In order to exclude subjects within the population who were 
genetically related, identity-by-descent estimates were performed for all possible 
combinations by PLINK vl.07 (see URLs). As a result, 2 XFS/XFG patients and 21 
controls who were assumed to be genetically related were excluded from the GWAS 
population. Ultimately, the population ended up with 201 XFS/XFG patients and 697 
controls (Figure 1 A and Table 1 ). Population stratification was examined by principal 
component analysis using EIGENSTRAT software v3.0 (see URLs). The four 
HapMap populations (CEU, YRI, JPT, and CHB) were simultaneously applied into 
EIGENSTRAT as the reference. The generated cluster plots indicated that the XFG 
and control population were genetically clustered within the JPT population, and 
there was no outlier sample (Supplementary Material, Figure S8). In order to extract 
high-quality autosomal SNP genotype data, the following QC filters were applied: (i) 
call rate per SNPs in each case and control samples ^95%, respectively, (ii) minor 
allele frequency (MAF) in case and control samples ^1%, and (iii) Hardy- Weinberg 
equilibrium (HWE) in control samples with a P-value ^0.001. Consequently, the 
remaining 652,792 SNPs were analyzed for the XFS/XFG GWAS (Figure lA). 

Replication study. For the replication study, an independent population of 384 
samples consisting of 121 XFS/XFG patients and 263 controls was prepared (Table 1). 
Those samples were extensively checked in order to remove any familial 
relationships. In total, 41 SNPs were applied to the replication by referencing the 
GWAS results as follows: 34 SNPs that passed the Bonferroni's correction in our 
GWAS and 7 SNPs which were selected to be useful for confirming the LD region of 
15q24.1 locus. In addition, 2 cSNPs (rsl048661 and rs3825942) were assessed that 
have been reported by Thorleifsson et al.^, yet are not designed in SNP 6.0. 
Furthermore, 4 SNPs associated with POAG in our previous study"*" were included in 
this replication assay. In total, 47 SNPs were divided into each of the two groups of 30 
and 17 SNPs per assay plate, respectively, according to the manufacturer's 
recommendation, and genotyped simultaneously for 384 samples in each group by 
use of the MassARRAY iPLEX Genotyping system (Sequenom, Inc., San Diego, 
California, USA). AU of the resulting call rates of these SNPs were over 90%. 

Meta-analysis. A meta-analysis was performed for 41 SNPs, which were genotyped in 
both the GWAS and replication study, by use of the Mantel-Haenszel test^". Among 
these SNPs, rs2165241 and rs4243042 were found to be monomorphic for the XFS/ 
XFG subjects used in the replication study (i.e., the risk allele frequency was 100%). 
Therefore, those 2 SNPs were excluded from the Mantel-Haenszel test. Genotype 
imputation for GWAS was performed using MACH ver. 1.0.16 software (see URLs). 
The genotype data of JPT and CHB in the HapMap Phase 2^^ (release #24) were 
applied as the reference for haplotyping. After the QC was applied for the imputation 
results (MAF 2:0.01 and Rsq >0.7), association analysis of the SNPs for XFS/XFG was 
performed by use of inach2dat ver. 1.1.9 software (Windows version). Conditional 
analysis was assessed for 34 genome-wide significant SNPs using the combined 
samples from the GWAS and the replication study (i.e., 322 XFS/XFG patients and 
960 controls). Haplotype analysis was examined with the samples from the 
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replication study in order to assess the reported cSNPs in LOXLl. Both conditional 
and haplotype analyses were performed by use of the PLINK vl.07 whole genome 
association analysis toolset. 

Data management and statistical analysis. In order to manage and analyze all of the 
genotype data, our in-house Genoika Server System'**'^^'^^ was used. The Genoika 
Server System comes with PLINK, R program with some packages, EIGENSTRAT, 
and Haploview (see URLs for the details of each software) built in, and all of the 
analyses were performed by use of this system. In addition, Microsoft Office Excel 
2003 (Microsoft Corporation, Redmond, Washington, USA) was used for preparing 
the data sets and statistical analysis. The frequency of alleles in the case and control 
samples was compared by use of the basic allele test. The odds ratio (OR) and the 
upper and lower limit of the 95% confidence interval (CI) of each SNP were calculated 
for the allele possessing a higher frequency m the case samples than in the control 
samples. The HWE was evaluated by use of the chi-square test. A regional plot in 
Figure 2 was drawn by use of LocusZoom ver. 1.1 (see URLs) genetic analysis 
software. Q-Q plots were generated by ranking the observed values from minimum to 
maximum and plotting them against their expected chi-square values using the 
"snpStats" package ver. 1.10.0 in the R program (see URLs). 

URLs. PLINK software vl.07, http://pngu.mgh.harvard.edu/"-purceLl/plink/; 
EIGENSTRAT software v3.0, http://genetics.med.harvard.edu/reich/Reich_Lab/ 
Welcome.html; MACH ver. 1.0.16 and mach2dat ver. 1.1.9 softwares, http://www. 
sph.umich.edu/csg/abecasis/MACH/index.html; Haploview v4.2, http://www. 
broadinstitute.org/scientific-community/science/programs/medical-and- 
population-genetics/haploview/haploview; LocusZoom version 1.1, http://csg.sph. 
umich.edu/locuszoom/. R programs, http://www.r-project.org/; SIFT dbSNP DB 
(buUd 132), http://sift.jcvi.org/; PolyPhen-2 version 2.2.2 r398, http://genetics.bwh. 
harvard.edu/pph2/. 
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